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ABSTRACT: We study the physical mechanism of how an external magnetic field influences the 
QCD quark condensate. Two competing mechanisms are identified, both relying on the interaction 
between the magnetic field and the low quark modes. While the coupling to valence quarks enhances 
the condensate, the interaction with sea quarks suppresses it in the transition region. The latter 
"sea effect" acts by ordering the Polyakov loop and, thereby, reduces the number of small Dirac 
eigenmodes and the condensate. It is most effective around the transition temperature, where the 
Polyakov loop effective potential is flat and a small correction to it by the magnetic field can have 
a significant effect. Around the critical temperature, the sea suppression overwhelms the valence 
enhancement, resulting in a net suppression of the condensate, named inverse magnetic catalysis. 
We support this physical picture by lattice simulations including continuum extrapolated results 
on the Polyakov loop as a function of temperature and magnetic field. We argue that taking into 
account the increase in the Polyakov loop and its interaction with the low-lying modes is essential 
to obtain the full physical picture, and should be incorporated in effective models for the description 
of QCD in magnetic fields in the transition region. 
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1 Introduction 

It is by now well established that strongly interacting matter has at least two distinct forms of 
existence: the hadronic and the quark-gluon plasma phase. They are separated by a crossover [1], 
occurring at a transition temperature T c of about 150 MeV [2-4]. Understanding the thermodynam- 
ics of strongly interacting matter around this transition is essential for the description of systems 
including the early universe, the interior of neutron stars or heavy-ion experiments. In heavy-ion 
collisions, for instance, the colliding beams can create extremely high magnetic fields at the colli- 
sion center, which can influence the thermodynamic properties of the system [5-8]. Similarly, the 
external magnetic field is expected to play an important role for very dense neutron stars (see, e.g., 
ref. [9]) and for the evolution of the early universe (see, e.g., ref. [10]). 

A particularly pronounced effect of the external field B is its influence on chiral symmetry 
breaking. At zero temperature, both low energy effective theories (see, e.g., refs. [11-15], and 
ref. [16] for a review) and early lattice simulations predicted the magnetic field to enhance the 
quark condensate. Although these first lattice simulations used the quenched approximation [17] and 
larger-than-physical quark masses [18, 19], the effect was later also confirmed by lattice simulations 
at physical quark masses in the continuum limit [20-22]. We can thus safely conclude that at 
zero temperature, the external magnetic field enhances the condensate. This phenomenon is called 
magnetic catalysis. We remark that the zero temperature magnetic catalysis of the condensate can 
also be derived from the positivity of the scalar QED /3-function, thereby showing its universal 
nature [23]. 

At finite temperature, and especially around T c - in the temperature range, which is most 
important for heavy ion collisons - the situation is not so clear. Here, most of the low energy 
models also predict that magnetic catalysis takes place, and, furthermore, T c is shifted to higher 
temperatures (for exceptions see, e.g., refs. [24, 25], and ref. [26] for a bosonic system with charged 
pion condensates). This behavior was also observed for larger-than-physical quark masses on coarse 
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lattices [18, 27]. In contrast, simulations on finer grids and with physical light quark masses indicated 
that around the crossover temperature exactly the opposite happens: the magnetic field suppresses 
the quark condensate and shifts the crossover temperature downwards [20-22]. The discrepancy 
between the different lattice results suggests that the use of fine grids and light enough quarks [20] 
might be essential for correctly describing the effect. Moreover, the remarkable similarity between 
the reaction of gluonic and fermionic observables to the magnetic field [28] implies that the indirect 
effect of B on the gauge degrees of freedom - exerted through the electrically charged quarks - must 
be a crucial ingredient for the description of the transition region. 

In the present paper, we study the physical mechanism behind the suppression of the condensate 
by the magnetic field in the transition region - an effect we refer to as inverse magnetic catalysis 1 . 
We will show that the magnetic field influences chiral symmetry breaking through two different 
mechanisms: one enhances and the other suppresses the condensate. Let us first describe the 
mechanism that enhances the condensate and is also responsible for catalysis at zero temperature. 
According to the Banks-Casher relation [30], the quark condensate is proportional to the spectral 
density of the Dirac operator around zero. The Dirac operator explicitly contains the magnetic field 
and, thus, its spectral density in a fixed gauge background depends on B. As we will show, the 
magnetic field enhances the spectral density around zero and, therefore, also enhances the quark 
condensate. In the free case as well as most of the model calculations this enhancement is linked 
to the degeneracy of eigenvalues being proportional to the magnetic flux. This is purely a "valence" 
effect that can already be observed in the quenched approximation, where the back-reaction of the 
quarks on the gauge field is ignored. We can view this as the basic mechanism explaining magnetic 
catalysis. 

Besides the explicit dependence of the Dirac operator on B, the magnetic field also enters the 
quark action, and influences the probability of different gauge backgrounds in the path integral. 
We find that this dependence leads to the suppression of the condensate in the transition region. 
This is a "sea" effect, since its origin is the different sampling of the gauge fields due to the quark 
determinant. The distinction between valence and sea effects was initially introduced in ref. [19] 
with a slightly different terminology (the "sea" was called "dynamical"). 

Although the two mechanisms are of very different nature, we will show that both of them orig- 
inate from the reaction of the low-lying Dirac modes to the magnetic field. However, the magnetic 
field "in the sea" has exactly the opposite effect in the transition region as the "valence" mechanism. 
Thus, the relative strength of the two effects will determine whether the QCD condensate under- 
goes magnetic catalysis or, rather, inverse magnetic catalysis 2 . As it turns out, the most important 
control parameter that influences the competition of the valence and sea effects, is the quark mass 
m. In particular, the two effects are expected to have a drastically different m- dependence. On 
the one hand, the valence effect is already seen in the quenched approximation, and is expected to 
depend on m only mildly. On the other hand, the sea effect is completely absent in the quenched 
approximation - which technically corresponds to infinite quark masses - and is expected to increase 
as the quarks become lighter. This shows that the quark masses have to be tuned to their physical 

x We note that the expression 'inverse catalysis' is also used for the reduction of the transition temperature with 
growing B at nonzero chemical potential [29], but presumably has a completely different physical origin. 

2 We note that possible explanations of the lattice results [20-22] regarding inverse magnetic catalysis were given 
recently in Refs. [31, 32] relying on the dimensional reduction of the lowest Landau level dynamics for large magnetic 
fields, eB > Aq CD , plus screening and neutral meson effects. We remark, that the lattice data indicate that inverse 
magnetic catalysis is related to the transition taking place, since for even higher temperatures T>T C the conventional 
magnetic catalysis again prevails [22]. 
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values in order to have a proper description of magnetic (inverse) catalysis. 

Besides identifying the sea mechanism as the one responsible for inverse magnetic catalysis, 
it is instructive to look at how the magnetic field in the determinant changes the typical gauge 
configurations contributing to the path integral. Understanding this is clearly essential for the 
incorporation of the sea effect into effective models describing QCD. As we will show, one of the most 
striking influences of the magnetic field on the gauge configurations is to drive up the expectation 
value of the Polyakov loop 3 . This can be understood by noting that the Polyakov loop has a strong 
influence on the lowest part of the Dirac spectrum. Namely, the ordering of the Polyakov loop shifts 
up the small Dirac eigenvalues and thus suppresses the quark condensate. However, the Polyakov 
loop is also influenced by the effective action originating from the gauge action and the fermion 
action at zero magnetic field. The change in this effective action caused by the magnetic field is 
usually relatively small, and cannot be decisive. The only exception to this is around the transition 
temperature, where the Polyakov loop effective action is flat, and a small contribution coming from 
the magnetic field can have a large ordering effect on the Polyakov loop. This, in turn, suppresses 
low Dirac eigenmodes and the condensate. We conjecture that this is the physical mechanism behind 
inverse magnetic catalysis seen around the transition temperature. 

The rest of the paper is devoted to a demonstration of the physical picture we just described us- 
ing numerical lattice simulations. The results we present are partly based on large scale simulations 
already reported in refs. [20, 22], and partly on new simulations involving small lattices, where a 
full diagonalization of the Dirac operator was feasible. In each case, we use 1 + 1 + 1 flavors of stout 
smeared staggered quarks with physical quark masses and the tree-level improved Symanzik gauge 
action. Details of the simulation setup and the configurations can be found in refs. [20, 35, 36]. 

The paper is organized as follows. In the next section, we define the valence and sea contribution 
to the quark condensate, and measure them for several magnetic fields at temperatures below and 
close to the transition. In Sec. 3, we discuss the effects of the magnetic field on the Polyakov loop, 
and Sec. 4 contains our conclusions. In the appendices we calculate the interplay of magnetic field, 
Polyakov loop and mass in the free case, and discuss the renormalization properties of the sea and 
valence condensates and the Polyakov loop. 

2 Sea and valence quark effects 

As we already noted, the most important notion, that connects all different aspects of the problem, 
is the spectrum of the quark Dirac operator. The Dirac spectrum is affected by the magnetic field for 
two different reasons. On the one hand, the magnetic field explicitly appears in the Dirac operator, 
and directly influences its spectrum in any fixed gauge background. On the other hand, through 
the quark determinant, B also affects the probability measure of how the gauge field configurations 
are sampled in the path integral. Throughout this paper, we will call the first one the "valence" 
effect, and the second one the "sea" effect. In this section, we introduce our notation, and discuss 
these effects in detail. 

In physical terms, the most important measurable property of the Dirac spectrum is the spectral 
density p(\) around zero, which is proportional to the chiral condensate due to the Banks-Casher 
relation [30], 



lim lim p(A) • tt = lim lim ipifj 
A— >0 V— >-oo m— >0 V— >oo 



TdlogZ 
V dm 




3 That the Polyakov loop plays a role for the catalysis mechanism in the transition region has been advocated in 
Refs. [33, 34]. 
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Here the condensate is denned in terms of the partition function, which is given by the path integral 
over gauge configurations U as 

Z{B) = [ VU e~ S9 det($(B) + m), (2.2) 



where S g denotes the gauge action. For the sake of clarity, here we consider only one fermion 
flavor with charge q and mass m, and suppress factors of 1/4, which appear due to the rooting 
procedure for staggered quarks. The temperature and the three-volume are given as T = (aA^) -1 
and V = (aN s ) 3 with N s (Nt) the number of lattice sites in the spatial (temporal) direction, and a 
the lattice spacing. We remark that since the magnetic field couples only to the electric charges of 
quarks, it enters exclusively in the combination qB. 

Expanding the derivative in Eq. (2.1), the condensate is obtained as 

^(B) = _i_ fvUe- Ss det(lp(B) + m)Ti(Ip{B) + m)- 1 , (2.3) 
Z ( B ) J 

showing that the magnetic field indeed appears both in the determinant and in the operator itself. 
To separate these dependences, we define the valence and sea condensates as 



^ val (B) = J VU e- s * det(lp{0) + m)Tiilp{B) + m)-\ 

^ sea (B) = — [ VU e~ s o det(lp(B) +m)Tr(#(0) + m)" 1 . 
Z \ B ) J 



(2.4) 



We note that valence condensates can be used to define dressed Wilson loops [37] , which are directly 
related to the QCD string tension in the large mass limit. 

Any physically consistent theory has to have the same valence and sea fermion content. Thus, 
at first sight it does not seem possible to separate the valence and sea effects of the magnetic field 
in a well-defined theory. One can, nevertheless, exactly reproduce the condensates in Eq. (2.4) 
by alluding to techniques from partially quenched QCD [38-40] (and quenched disorder [41, 42]). 
Using commuting spin 1/2-fields (so-called ghost quarks), one can generate inverse Dirac determi- 
nants in the functional integral. With adjusted charges and masses, these inverse determinants can 
cancel 'unwanted' determinants in the path integral, arriving at the valence and sea condensates 
of Eq. (2.4). From a different point of view, one can also directly obtain the sea condensate in a 
theory with an electrically charged and a neutral fermion flavor, by looking at the condensate of the 
neutral fermion in the presence of the magnetic field. Since B appears in the determinant of the 
charged flavor, but not in the neutral propagator, this indeed isolates the sea effect. Even though 
in QCD all fermion species are electrically charged, on the technical level of the lattice theory, the 
valence and sea effects are naturally separated. We will use a similar argument in App. A to discuss 
the renormalization of the sea and valence condensates. 

To discuss the effect of the external magnetic field, we are interested in the change of the 
condensates due to a nonzero B. This change is given by the difference 

2tt7 — — 
AS (^) = T^2 [WW - ^(0)] • (2-5) 

This combination is particularly useful, because both additive and multiplicative divergences cancel 
in it. Based on the Gell-Mann-Oakes-Renner relation, the normalization is chosen such that AS is 
measured in units of the condensate at zero magnetic field and zero temperature [22]. We define 
A£ val and A£ sca in a similar manner from Eq. (2.4). At B = the three types of condensate 
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Figure 1. The valence (left panel) and sea (right panel) contributions to the up quark condensate as a 
function of the magnetic field, calculated at two different temperatures. The two temperatures are chosen to 
be well below and around T c . The different symbols correspond to three different lattice spacings (decreasing 
as N t grows). 



are obviously equal. Furthermore, for small magnetic fields (assuming analyticity in qB), the two 
contributions appear additively in the total condensate [19], 

AE(S) ~ A£ val (£) + AZ sca (B), (2.6) 

showing that this separation indeed makes sense, at least for small magnetic fields. In App. A we 
show that A£ val and AS sca are both properly renormalized. 

In Fig. 1, we show how the valence condensate AE val (left panels) and the sea condensate A£ sea 
(right panels) for the up quark depend on the magnetic field, at two different temperatures. The 
two temperatures were chosen to be well below and just below the transition temperature. Clearly, 
the magnetic field in the valence Dirac operator enhances the condensate at both temperatures. In 
contrast, the sea effect enhances the condensate only well below T c , whereas around T c it suppresses 
it. Eventually - around T = 160 MeV - the sea contribution becomes the dominant one, resulting 
in a decrease in the total condensate and thus inverse magnetic catalysis [22]. Fig. 1 contains 
results for three different lattice spacings, showing that the effect persists in the continuum limit 
as well. This is in sharp contrast to the findings of ref. [18], where an enhancement was found at 
all temperatures, for larger-than-physical quark masses. In Sees. 2.1 and 2.2, we discuss how the 
underlying mechanism responsible for the valence and the sea effects can be understood, based on 
eigenvalues of the Dirac operator at nonzero magnetic fields. 

2.1 The valence effect 

The valence effect can be easily understood by inspecting how the low part of the spectral density 
p(X) of the Dirac operator depends on the magnetic field in any gauge field background. To be 
specific, we use a set of gauge field backgrounds generated at zero magnetic field, but the qualitative 
picture is the same in any reasonable ensemble of gauge fields. 

In Fig. 2, we plot the spectral density of the Dirac operator for three different values of the 
(valence) magnetic field, as measured on Nt = 6 lattices, generated at B = and T = 142 MeV. 
This gauge field ensemble corresponds to the N = 6, T = 142 MeV data for the valence condensate 
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Figure 2. The spectral density of the Dirac operator around zero, computed at three different values of 
the magnetic field. In all three cases, the gauge configurations were generated without the magnetic field in 
the quark determinant. 



in Fig. 1. One can clearly see the increase of the spectral density and, thus, of the valence con- 
densate with the magnetic field. We remark that the same behavior is reproduced for any gauge 
background, independent of the temperature and magnetic field, which was used for the generation 
of the configuration. In other words, this means that the change in the valence condensate AS val 
is always positive. We note that a similar proliferation of low Dirac eigenmodes already occurs in 
the free theory, see the discussion in App. B. Moreover, a remarkable feature of the free spectrum 
on the lattice is that the eigenvalue pattern as a function of the magnetic field is similar to the 
so-called Hofstadter butterfly, the energy levels of Bloch electrons in a magnetic field [43] . 

2.2 The sea effect 

The sea effect arises, because the magnetic field in the quark determinant changes the relative 
weight of the gauge configurations, and is therefore equivalent to a reweighting in B. In general, 
reweighting is a technique that uses configurations generated at a given (starting) point of the 
parameter space, and assigns a new weight to each configuration, in a fashion that the resulting 
ensemble describes the system at a new (target) point of the parameter space. Thus, the expectation 
value of an arbitrary observable O at the target point is obtained in terms of measurements on the 
configurations generated at the starting point. Here, we will consider the B = system as the 
starting point, and the B > system as the target ensemble. For this case, the difference of weights 
equals the ratio of quark determinants at B and at B = 0, and the exact rewriting of the expectation 
value at B reads 



(O) 



1 



B 



m 

Z(B) Z(0) 
e -AS f (B) 



VUe 



, det(Ip(B) + m) 

det (00 + m) Z/ / + 

V V ' ' det(0 O +m 



O 



e -AS f {B)\ 
I 



(2.7) 



where the subscript of the expectation value indicates the value of the magnetic field, at which the 
ensemble is generated. Here, we defined the change in the fermionic action due to the magnetic 
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field, 

- AS f (B) = log det(J0(B) + m) - log det(0(O) + m). (2.8) 

Let us now apply this machinery for the case of the condensates. The valence condensate of 
Eq. (2.4) equals the B = expectation value, 

^ val (^) = <Tr(0( J B)+m)- 1 > o , (2.9) 

whereas the full condensate is obtained in the reweighting picture as 

W>(B) = (e~^ s ^ Tr(]p(B) + m)" 1 ) / (e~ A5/(B) ) Q • (2.10) 

Now we are in the position to discuss the sea effect, i.e. the difference between the full and valence 
condensates. The full condensate is written in Eq. (2.10) as the product of the trace of the prop- 
agator at nonzero B and the factor exp(— ASf(B)), both of which depend on the configuration. 
Therefore, we have to study the correlation between the fermion action difference and the value 
of the condensate at nonzero B. To this end, in Fig. 3 we show a scatter plot of the down quark 
condensate versus ASf(B) for an ensemble of 150 10 3 x 4 lattices generated at B = 0, and at inverse 
gauge coupling (5 = 3.35, corresponding to a temperature around T c . On this gauge ensemble, we 
computed the condensate by inserting eB ~ 0.5 GeV 2 in the Dirac operator. Both the action and 
the condensate were calculated in the spectral representation, by explicitly diagonalizing the quark 
Dirac operator. 
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Figure 3. Scatter plot of the down condensate versus the change in the fermionic action due to the 
magnetic field, Eq. (2.8), at eB « 0.5 GeV 2 . Each point represents a 10 3 x 4 gauge configuration. 

In Fig. 3, each point represents a gauge configuration. The ordinate is the condensate and the 
abscissa is the action difference on that particular gauge configuration. Since the configurations were 
generated at B = 0, a simple averaging of the ordinates in Fig. 3 would give the "valence" condensate, 
in accordance with Eq. (2.9). To obtain the full physical condensate, including the valence and sea 
effects, one needs to weight each ordinate with the Boltzmann factor exp(— ASf). This means that 
the contribution of points corresponding to larger abscissas are exponentially suppressed. A clear 
tendency can be seen in the figure, showing that larger values of the condensate generally correspond 
to larger actions, and are therefore suppressed by the quark determinant. This is the "sea" effect. 

The reweighting in the magnetic field can substantially influence the condensate, because the 
reweighting factors ASf(B) fluctuate sufficiently strongly. In Fig. 3, for example, a distance of 7 
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on the abscissa amounts to a change in the path integral weight by three orders of magnitude. It 
is known that the light quark condensate is dominated by low Dirac modes. We now show that 
the fluctuations of the change in the quark action also have their origin in the low spectral region 

of the Dirac operator. In Fig. 4, we plot how the variance o(x) = (x 2 ) — (x) 2 of the change in 
the fermion action x = ASf(B) builds up starting from the low-end of the Dirac spectrum. We 
approximate the determinants by the product of the lowest i eigenvalues, and plot the variance of 
this "approximated" ASf(B) as a function of i, on the same gauge ensemble that was used above. 
Looking at the figure, the dominance of the low Dirac modes in this quantity is obvious. 

We emphasize that these reweighting factors encode all information 4 about how the system 
evolves and how different observables change as the magnetic field grows. The magnetic field can 
have an influence on an observable only if two conditions are met: the reweighting factors have a 
non-zero variance and they are correlated with the given physical quantity. Since the variance of 
the reweighting factors comes mostly from the low Dirac modes, they are expected to be correlated 
with any other quantity that is sensitive to the low Dirac eigenvalues. In this context, the sea effect 
of the magnetic field on the condensate is naturally expected. 




Figure 4. The variation of the change in the effective action ASf (at eB ss 0.5 GeV 2 ), calculated using the 
lowest i eigenmodes, as a function of i. The total number of eigenvalues on this lattice is 6000. 



3 Connection to the Polyakov loop 

We have seen that the quark determinant tends to suppress gauge configurations with larger values 
of the quark condensate. To understand the reason behind this, it is desirable to identify those gauge 
field degrees of freedom that play the most important role in this mechanism. Such a description 
also facilitates the incorporation of this effect into low-energy effective models. 

Around the transition temperature, the quantity depending most sensitively on the control 
parameters is the order parameter. In the case of QCD, the approximate gauge field order parameter 

4 On the technical level, due to finite statistics, the applicability of reweighting is limited to starting and target 
systems whose important configurations are not too different. This potential overlap problem is not relevant here, 
since we use reweighting to unravel a mechanism and the results can always be checked by direct simulations of the 
target system (nonzero B). 
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around T c is the (averaged real part of the traced) Polyakov loop, 

1 / \ 
P=-/^ReTr ]J C/ 4 (x,t)\. (3.1) 

\ x t=0 / 

In the valence/sea language the Polyakov loop is a purely "sea" observable, since the expression to 
be averaged over does not depend explicitly on B (this applies to every gluonic observable). 

To see how the Polyakov loop reacts to the magnetic field, in Fig. 5, we show a scatter plot 
of P versus the change ASf(B) in the fermionic action due to the magnetic field. This plot is 
similar in spirit to Fig. 3, and is based on the same gauge configurations. The two plots also have 
a qualitatively similar appearance, but in the case of the Polyakov loop, the correlation is opposite. 
Here a simple average over the ordinates would give the B = Polyakov loop, whereas a weighted 
average yields the Polyakov loop at nonzero B. The fermionic action, and hence the magnetic field, 
suppresses small values and favors large values of the Polyakov loop. Below we will argue, that large 
Polyakov loops point to smaller condensates. 
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Figure 5. Scatter plot of the Polyakov loop versus the change in the fermionic action due to the magnetic 
field at eB w 0.5 GeV 2 . Each point represents a 10 3 x 4 gauge configuration. 

It is also not hard to understand why this effect is so pronounced around the transition tem- 
perature. Since there is a crossover at T c , the effective potential for the Polyakov loop is flat there. 
As a result, even a small contribution coming from the magnetic field in the quark determinant 
can significantly change its expectation value. In the following, we explicitly demonstrate this by 
calculating the dependence of P on the external magnetic field in the continuum limit. To extrap- 
olate the Polyakov loop to the limit a — > 0, we perform its proper renormalization, considering 
the scheme suggested in ref. [44], which we generalize to the case of nonzero magnetic fields. This 
renormalization is of the multiplicative form 

P r (a,T,B) = Z(a,T)-P(a,T,B), (3.2) 

and is discussed in detail in App. C. 

For the Polyakov loop calculation, we again employ the gauge configurations of ref. [20], gen- 
erated with physical quark masses at various values of the temperature and the magnetic flux 
Nfj ~ eB/T 2 . In any finite volume, this flux is quantized [45]. In order to determine the Polyakov 
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Figure 6. Left panel: the temperature- and magnetic flux-dependence of the renormalized Polyakov loop 
in the continuum limit. The solid lines represent curves of constant magnetic field {eB ~ N^T 2 values as in 
the right panel) . Right panel: the dependence of P r on the temperature around the crossover region. The 
different types of curves indicate lattice results obtained with different lattice spacings (different temporal 
extents) . The shaded areas show the continuum extrapolations together with their uncertainty. 

loop as a function of T, along a line of constant eB, an interpolation between the different fluxes 
Nb is necessary. We carry out this interpolation in a systematic manner, by fitting our data points 
for all temperatures, magnetic fluxes and lattice spacings altogether by a lattice spacing-dependent, 
two-dimensional spline function. A similar spline fit is described in ref. [46]. Due to the scaling 
properties of the action we use, the dependence on the lattice spacing is expected to be quadratic. 
We incorporated this in the fit by having two parameters on each node point as p\ +P2 • a 2 . Taking 
eB = const, slices of this two-dimensional surface at a certain a gives the Polyakov loop for that 
particular lattice spacing, while the a = surface corresponds to the continuum limit. 

In the left panel of Fig. 6, we plot the continuum extrapolated renormalized Polyakov loop P r 
as a function of the temperature and the magnetic flux. The solid lines upon the surface correspond 
to eB = 0, eB = 0.45 GeV 2 and eB = 0.75 GeV 2 slices. In the right panel of the figure, we 
show the temperature dependence of P r for these magnetic fields on the three lattice spacings, 
together with the continuum extrapolation. The shaded bands represent here the uncertainty of 
the continuum extrapolated P r . The results clearly show, that the Polyakov loop increases sharply 
with the magnetic field around T c , and that this feature persists in the continuum limit as well. As 
an empirical finding from that figure, inflection points of these curves are not very precisely defined, 
but the transition temperature from the renormalized Polyakov loops clearly decreases with the 
magnetic field. 

In the previous sections, we saw that low Dirac modes are the key to understanding inverse 
magnetic catalysis. To complete the picture, we would like to discuss one more point, namely 
the relationship between the low Dirac modes and the Polyakov loop. It is well-known that light 
dynamical fermions break the Z(3) center symmetry of the quenched theory, by forcing the system 
into the real Polyakov loop sector. This can be most easily understood starting from a free field 
picture. If the gauge field background is trivial - apart from a spatially constant Polyakov loop - the 
lowest Dirac eigenmodes are constant in space, and change smoothly in the time direction to fulfill 
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the boundary condition. The boundary condition and the phase of the Polyakov loop then combine 
to give the effective boundary condition (twist). The lowest fermion eigenvalue A m i n - a generalized 
Matsubara frequency - is proportional to this temporal twist (see discussion around Eq. (B.2)). For 
the usual antiperiodic boundary conditions, the twist is maximized when the Polyakov loop is in 
the real sector. Therefore, A m i n is also maximized for real Polyakov loop. Since the difference in the 
quark action between the Polyakov loop sectors is dominated by the small Dirac eigenvalues [47], 
the real Polyakov loop sector, the one with the largest A m i n , will result in the largest determinant. 

It turns out that a qualitatively similar picture applies also to the interacting theory above the 
transition temperature: a larger temporal twist pushes the lowest modes higher up in the spectrum, 
away from zero [48-50]. As a result, the fermion determinant prefers the real Polyakov loop sector 
in the SU(3) case. This is how the Z(3) symmetry is broken by dynamical fermions. If the Dirac 
operator has a tendency of having more small modes, the ordering of the Polyakov loop can shift 
more small eigenvalues up, and the dependence of the determinant on the Polyakov loop is stronger. 
This results in a stronger symmetry breaking effect and a stronger ordering of the Polyakov loop. 
This is exactly what happens in the presence of the magnetic field. 

We stress that an overall shift in the fermion action by the magnetic field would not produce 
any physical effect. The important point is that this shift fluctuates, and its fluctuations correlate 
with physical quantities (Polyakov loop, condensate). This is what produces the reweighting effect 
that we discussed. We have demonstrated this by virtue of scatter plots in Figs. 3 and 5. For a 
qualitative understanding, let us extend the Matsubara picture by magnetic fields. As the calculation 
in App. B shows, the free energies at nonzero B favor larger, i.e. deconfined Polyakov loops. Here, 
we focus on the main effect for large B, where only the lowest Landau level is occupied. Its 
eigenvalue is independent 5 of B, but its degeneracy still increases proportional to the magnetic flux. 
The lowest mode A m j n of a given configuration thus acts as a "handle" for the magnetic field to 
change the fermionic determinant, to a first approximation proportional to (Xmin) K suppressing 
configurations with small eigenvalues even further. At the same time, configurations with many 
small eigenvalues are correlated to small Polyakov loops and large condensates, which are therefore 
suppressed as well. This is what is observed in full QCD around the transition. We repeat that 
the validity of this picture depends on low mode dominance and the sensitivity of eigenmodes and 
condensates to Polyakov loops. Both are expected to hold for temperatures above the transition. 
It would be interesting to study in detail how the structure of the low quark modes change due 
to the magnetic field. An exciting clue in this direction might be the presence of independent low 
modes at high temperatures, localized to Polyakov loop islands [51]. The magnetic field might have 
a more coherent effect on these localized modes than on the bulk modes, following random matrix 
statistics. However, a detailed study of this mechanism is beyond the scope of the present paper. 

We note that previously in the PNJL model [52] and also in lattice simulations with larger-than- 
physical quark masses [18] the Polyakov loop was observed to decrease with the magnetic field, and 
the transition temperature was shifted up by B. We conjecture that this is the reason why in those 
cases inverse magnetic catalysis was not observed. Inverse magnetic catalysis depends crucially on 
the interaction between the quark determinant and the Polyakov loop, which happens through the 
lowest part of the Dirac spectrum. This delicate effect has to be properly accounted for, to get the 
full physical picture. Since the Z(3) symmetry breaking of the quark determinant depends strongly 
on the quark mass, lighter quarks are expected to enhance inverse magnetic catalysis. Let us finally 
remark that the latest results of the authors of ref. [18] also indicate the sea condensate to decrease 

5 for spin 1/2 particles without anomalous magnetic moment 
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with B for higher temperatures T ~ 220 MeV, where the Polyakov loop is observed to increase with 
growing B [53], thus exhibiting the same mechanism as described above. 

4 Conclusions 

Using lattice simulations, we have analyzed how an external magnetic field affects the QCD vacuum 
at finite temperatures. Let us summarize: when the magnetic field is switched on, the quark 
determinant suppresses gauge configurations, on which the Dirac operator with the magnetic field 
has many small eigenvalues. Around and above T c , this suppression can happen by ordering the 
Polyakov loop. This ordering effect is particularly efficient around T c , where the Polyakov loop 
effective potential is flat, and the magnetic field in the determinant has a significant ordering effect 
on the Polyakov loop. An ordered Polyakov loop implies that the lowest Dirac eigenmodes become 
similar to Matsubara modes, and are shifted up in the spectrum. For small quark mass, the quark 
condensate gets the largest contribution from the low-end of the Dirac spectrum. Therefore, fewer 
low eigenmodes also imply a suppression of the quark condensate. This suppression of the condensate 
that we called "sea effect", competes with the enhancement of the condensate due to the magnetic 
field in the operator ("valence effect"). Due to the above described reasons, the sea suppression is 
particularly efficient around T c , and it actually overwhelms the valence enhancement, resulting in 
inverse magnetic catalysis. The unifying theme in our discussion is the lowest part of the Dirac 
spectrum. The response of low Dirac modes to the magnetic field can explain the valence and sea 
effects, as well as the ordering of the Polyakov loop. 

The low mode dominance, which we have used in our arguments, and illustrated in Fig. 4, 
strongly relies on the smallness of the quark masses. One might speculate that the effect of inverse 
magnetic catalysis becomes even more pronounced closer to the chiral limit. In any case, we believe 
to have identified an important mechanism at work in the transition region, already for magnetic 
fields slightly above the QCD scale (in contrast to asymptotically large ones). It would be useful to 
incorporate such effects back into low-energy effective models for QCD. Our results indicate that 
a proper account of the interaction between the Polyakov loop and the Dirac determinant and its 
influence on the condensate is needed for that. 
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A Renormalization in the valence and sea sectors 

In this appendix, we show that the valence and sea condensates AE val and AS sea are properly 
renormalized, following the argumentation of refs. [54, 55]. Let us write down again the partition 



- 12 - 



function, as in Eq. (2.2), for a single quark flavor with charge q±, 

Z = Z qiB = j VUe- s ^et{U){ qi B) + m), (A.l) 
from which the condensate is derived to be 

MB) = £^ff^ = I WCftS) + m)-X B , (A.2) 

where the subscript indicates that the expectation value is evaluated with respect to the partition 
function Z qi B- The additive divergences of the condensate cancel in the change A with respect 
to the magnetic field. To eliminate the multiplicative divergence, a further multiplication by the 
quark mass is necessary, as in Eq. (2.5). Altogether, the combination AS ~ mAipip is ultraviolet 
finite. On the other hand, in the ensemble described by Z qi B, the change in the valence and sea 
condensates of Eq. (2.4) read 



A£ val ~ mA# ral = ?T"» [(Tr(^(9i5) + m)- 1 ) - (Tr(#(0) + m)- 1 )J , 



V 



TV($(0) + m)-\ B - <Tr(0(O) + m)~ 



(A.3) 



Our objective is to show that these combinations are also properly renormalized. 

To this end, let us consider an additional quark flavor, with mass m 2) charge q 2 , and multiplicity 
e. Even though the multiplicity is in practice restricted to take integer values, in principle one can 
also define a theory with a real multiplicity e S R, described by the partition function, 

Z qi B + e Q2 B = J VUe- s * det (J0(qiB) + m) det {${q 2 B) + m 2 ) e . (A.4) 

The condensate of the additional quark flavor is correspondingly given by 

$ MB) = TdJo g Z^B = T £ {Tr{]p{q2B) + ^-^^ > (A . 5) 

where this time the expectation value is evaluated in the ensemble with both the original and 
the additional flavor taken into account. Let us now expand Eq. (A. 5) around e = 0. Since 
det e = 1 + 0(e), we get the expectation value in the Z qi B ensemble plus higher order corrections, 

= £ • e ■ (Tr(J0(q 2 B) + m^ 1 }^ + 0(e 2 ). (A.6) 

Again, the combination m2A'0 e '0 e is completely renormalized for any value of e. Thus, in its 
expansion with respect to e, all terms must also separately be ultraviolet finite. In particular, the 
0(e) term at mi = m 



T 

-■m 



Tr{U)(q 2 B) + m)- 1 )^ - (Tr(#(0) + m)- 1 >J (A.7) 



is properly renormalized, for any values of the charges q\ and q 2 - 

Upon substituting q\ = and renaming q 2 = qi, Eq. (A.7) exactly equals the valence condensate 
of Eq. (A.3). On the other hand, substituting q 2 = gives the sea condensate. Altogether, this 
shows that the combinations AS val and AS sea , as defined after Eq. (2.5), are properly renormalized. 
The above argumentation relies essentially on the fact that the magnetic field enters only in the 
combination qB. Therefore, varying the charges effectively allows for varying the magnetic field. 
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B Free case calculation 



We consider the free energy of (one flavor of) massive quarks in the continuum at finite temperature 
and study the interplay of constant Polyakov loops, mass and magnetic field (and chemical potential 
as a cross-check). 

Let the gauge background be an abelian field generating a constant magnetic field (0, 0, B), e.g. 
a x = 0, a y = Bx, plus a constant SU(3) gauge field Aq generating a diagonal and real Polyakov 
loop, 

( e 2niip \ 
e~ 2 ^ J , l TrL = l (l + 2cos(2vrv?)). (B.l) 

Obviously, ip = 1/3 amounts to a traceless Polyakov loop, which we will take as a signal for the 
confined phase, while for ip — > the Polyakov loop becomes trivial, L — >• 1, which stands for the 
deconfined phase. 

In the computation of the eigenvalues of the Dirac operator ifi = Ta-Du with covariant derivative 
Dn = dn + iA^ + iqdfj, at finite temperature T, the Polyakov loop background can be incorporated in 
the quark boundary conditions in each color sector as t/j a (xo + l/T, x) = — exp(27rix a )V'a(2 ; 0) with 
the three phases Xa £ {</>> 0}. This modifies the Matsubara frequencies to (it + 2-nXa + 2irk)T 
with integer k. 

The eigenvalues of the squared Dirac operator in each color sector read 

- U) 2 a {B) -> X 2 a (B) = (2tt(1/2 + X a + k)Tf + p 2 z + |g£|(2n + 1) + gS S; (B.2) 

with momenta p z E R, Landau levels n = 0, 1, . . . and s = ±1 being twice the spin. The degeneracy 
of these eigenvalues is proportional to the magnetic flux B ■ A with A the area of the system 
perpendicular to the magnetic field. Together with the z-extension L z in the p z integration, this 
yields the three-volume V = AL Z in 

A 2 V ; fc=-oo J n=0 s=±l 

The factor 2 comes from the degeneracy of particles and antiparticles at given spin (both together 
provide the four eigenvalues of the four-dimensional gamma matrices). 
For the partition function we make use of chiral symmetry to write 

3 1 
logZ(fl) = X>gdet(# a (5) + m) = - ^ ^ log(A^(5) + m 2 ). (B.4) 

a=l a A 2 

The proper time/zeta function regularization formalism can be introduced via the Mellin transform 

d 



logy = 

oa 



y = - bt 



>o r ( a ) 



1 f°° 

— / dss^eMsy), (B.5) 



(valid for Re y, Re a > 0). Through interchanging the order of limits, integrals etc. analytic 
continuation is performed and (some of) the singularities are removed. This technique has been 
used in the classical papers [56, 57], for a review see e.g. [58]. 

Using Eqs. (B.3-B.4), the free energy density / = — TlogZ/V for a given background becomes 

f(tp; B) = ^J2j o ds s ~ 2 exp(-m 2 s) coth(qBs) 6 3 (( Xa + 1/2)tt, e"5^) , (B.6) 
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where we have assumed positive qB for simplicity and performed the Gaussian integration over p z , 
the derivative wrt. a and the sums over n, s and k. The latter introduced an elliptic theta function 
9 3 (u, q) = l + 2 Zn=i <f 2 cos(2nn). 

The s-integral is known to be divergent for small s (as coth(qBs) — > l/(qBs) and . . . , exp( — )) 
1). Since s has dimension length squared, this is a UV-divergence. After subtracting the B = 
part it is weakened to be logarithmic and can be absorbed in the renormalization of the electric 
charge [57]. 

Here we follow a different route and consider the difference of the free energy of a given back- 
ground (p and the deconfining background (p = 0, which means 

A v f(B) = f(tp; B) - /(0; B) (B.7) 
= ~fJ o dss~ 2 exp(-m 2 s)coth(qBs) 6> 3 ((v? + l/2)vr, e~3^) - #3 (n/2, e~3^) 

where we used the ip — > —ip symmetry. This integral is finite because the nonperturbative behavior 
exp(— #/s) of the theta function difference for small s now kills all power-like singularities. 

Note, that the integrand in (B.7) is positive definite and so is A^f. Hence the quark determi- 
nant always favors ip = 0, i.e. deconfining Polyakov loops. Furthermore, the integrand increases 
monotonically as m is reduced, such that this effect becomes more pronounced for light quarks. 
This is consistent with the finding from 'lattice experiments' that the onset of the deconfined phase 
- the deconfinement transition - occurs at lower temperatures, when the quarks become lighter. 

As another cross-check we consider the effect of an imaginary chemical potential [i = if]. It 
can also be incorporated in the boundary condition (for all color sectors in the same way) and the 
corresponding change of free energies follows from Eq. (B.7) in the B — > limit is 

1 f°° r / _ 1 \ 

= 2vr2 / dS S ~ 3 exp (~ m2,s ) l 03 [if + 1 / 2 ~ ? ?) 7r ' e (B.8) 

J 

+ 3 ((-¥> + 1/2 -77)7r,e~5^) -2 3 ((1/2 - rfiir, e"i^ 

Again, the integrand is positive definite, and the effect of the quark mass is the same. Increasing 77 
we have observed A„/(r/) to decrease, although the integrand does not change monotonically with 
rf. Thus an imaginary chemical potential r\ favors <p > 0, i.e. confining Polyakov loops, consistent 
with the finding that the transition temperature grows with it. 

After these checks, we analyze the effect of an external magnetic field in the same fashion. For 
this the dependence of the integrand in Eq. (B.7) on qB is decisive: the function qB ■ coth(qBs) 
increases with qB for all s. Hence a magnetic field favors deconfined Polyakov loops and tends 
to lower the transition temperature. This is certainly a sea effect as it comes from the quark 
determinant in the effective action, and it gets washed out by large masses m. This effect is already 
present in the lowest Landau level (n = 0, s = —1 in Eq. (B.2)). The result of this projection is 
to replace coth(qBs) — > 1 in Eq. (B.7). The integrand still increases proportional to qB, from the 
degeneracy of the modes. 

It is interesting to note that these findings are inverted when modifying either the statistics, 
the boundary condition or the spin of the quarks: for commuting c-number quarks the determinant 



6 The first derivatives of A v f(n) with respect to rj and ip vanish at 77 = ip — 0, i.e. to linear order in the imaginary 
chemical potential and around the deconfined phase. The second derivatives with respect to rj and ip at T) = (p = 
give T 2 with a negative coefficient, i.e. the decrease mentioned. 
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has to be replaced by its inverse giving a minus sign to the free energy. For quarks periodic in the 
time-like direction one has to remove tt/2 in both 63 s and the integrand becomes negative definite. 
For spinless quarks qB ■ coth(qBs) gets replaced by qB/ smh(qBs) which decreases with qB. 

We conclude this section by a brief discussion on the validity of these considerations. First of 
all, the approximation of QCD by free quarks is expected to work at high temperatures thanks to 
asymptotic freedom. The free energy for instance, approaches the Stefan-Boltzmann result and the 
spectral gap has been found to be described by the lowest eigenvalue in the free case [51, 59] (with 
background Polyakov loop as we do here). 

The (de)confinement order parameter is the spatially averaged Polyakov loop. It is this quantity 
that changes from zero to non-zero (small to large in a crossover) as the temperature rises. The local 
Polyakov loops strongly vary, as can be seen from the decaying Polyakov loop correlator (related 
to one definition of confinement) and the distribution of local Polyakov loops (which in the low 
temperature phase is known to be the Haar measure). Such a local resolution of the Polyakov 
loop is lost in our approach as we compare spatially constant backgrounds. Thus, we cannot speak 
about transitions between disordered and ordered Polyakov loops. Nevertheless, the finding that 
free quarks favor constant trivial Polyakov loop backgrounds in comparison to constant nontrivial 
ones can be interpreted as a signal for the suppression of the confined phase by quarks in realistic 
QCD. This effect has the right tendency under a change of the control parameters quark mass and 
(imaginary) chemical potential. Moreover, as we have demonstrated, for external magnetic fields it 
correctly describes the increase of the Polyakov loops and the decrease of the transition temperature 
as seen in full QCD. 

C Renormalization of the Polyakov loop 

In this appendix we discuss the renormalization of the Polyakov loop P. First of all, we note that P 
is related to the free energy of an infinitely separated, static quark-antiquark pair F qq = —2TlogP, 
which contains additive divergences in the cutoff. These additive divergences, therefore, appear as 
multiplicative divergences in P. A possible scheme for the corresponding renormalization is to set 
the renormalized free energy to a fixed value F±, at B = and a fixed but arbitrary temperature 
T± > T c . Since the magnetic field introduces no new divergences in F qq (see below), this subtraction 
is unchanged at B > 0, making the combination 

F qq , r (T, B) = F qq (T, B) - F qq (T„ B = 0) + F* (C.l) 

finite in the continuum limit. 

To see that the divergences of F qq are indeed ^-independent, let us consider the total free energy 
F of the system, in the presence of an external magnetic field. F includes both the contribution 
of matter, and that of the field itself, B 2 /2. It has been known for long [57] that the only in- 
dependent divergence, which resides in F(B), is eliminated via the simultaneous renormalization 
of the electric charge q, and that of the pure magnetic term B 2 /2 (see also refs. [23, 58]), which 
leaves the combination qB invariant. On the other hand, the free energy F qq of the static color 
charges only couples to the magnetic field via sea quark loops, in the form qB, and contains no pure 
magnetic term. Then, it also cannot contain any independent divergence. 

One can also explicitly show the absence of independent divergences in F qq using perturbation 
theory. Namely, the leading order perturbative term in F qq is given by the exchange of a gluon. 
The latter creates a virtual (sea) quark loop, to which B couples. This diagram is the same as 
the quark loop with two photon lines, up to constant group theoretical factors. Using the exact 
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quark propagator in the magnetic field [57], this photon self-energy has been explicitly calculated 
at asymptotically large magnetic fields to be finite [60] . Diagrams corresponding to higher orders in 
perturbation theory contain the same diagram as subdiagram, or sea quark loops with even more 
gluon lines attached, which (since they contain a higher number of quark propagators) are also 
finite. Altogether this shows again that the renormalization of Fq q - and, accordingly, that of P - 
is independent of B. 

We can now translate the renormalization condition (C.l) for Fq q to the renormalization of 
the Polyakov loop. On a lattice with spacing a, the corresponding renormalization for P = 
exp(—Fq q /2T) amounts to 



where P+ = exp(— F*/2T+). We choose T* = 162 MeV, where the bare Polyakov loops are already 
significantly nonzero, and we set F± = (P* = 1). Different values for these parameters correspond 
to different renormalization schemes, connected by finite renormalizations. 
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